How unique is the expected stress-energy tensor of a massive scalar field? 
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We show that the set of ambiguities in the renormalized expected stress-energy tensor allowed by 
the Wald axioms is much larger for a massive scalar field (an infinite number of free parameters) 
than for a massless scalar field (two free parameters). We also use the closed-time-path effective 
action formalism of Schwinger to calculate the expected value of the stress-energy tensor in the 
incoming vacuum state, for a massive scalar field, on any spacetime which is a linear perturbation 
off Minkowski spacetime. This result generalizes an earlier result of Horowitz and also Jordan in 
the massless case, and can be used as a testbed for comparing different calculational methods. 
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I. INTRODUCTION AND SUMMARY 

A. Background and Motivation 

In semiclassical gravity, a classical metric is coupled 
to quantum fields according to the semiclassical Einstein 
equation 



Gab — 8nG (T a b 



(1.1) 



where G a b is the Einstein tensor and G is Newton's grav- 
itational constant. This equation is usually postulated 
rather than derived as there is no complete theory of 
quantum gravity from which it could be derived, al- 
though several formal derivations have been given |Q. 
There are several well-known difficulties associated with 
the semiclassical theory. First, there are difficulties as- 
sociated with the existence of unphysical, exponentially 
growing "runaway" solutions of Eq. ( p-.l| ) , which have not 
yet been completely resolved The second difficulty, 
which is the subject of this paper, is the non-uniqueness 
of the expected stress-energy tensor on the right hand 



side of Eq. (1.1 



For a scalar field, several methods have been suggested 
for calculating the expected stress energy tensor. These 
include (i) the "point splitting" algorithm ||, (ii) the 
deWitt-Schwinger expansion method 0, and (iii) the 
closed-time-path or in-in effective action method HQ. 
There is no general agreement as to which method is 
correct. For example, it is claimed in Ref. || that the 
deWitt-Schwinger method is invalid for a massive scalar 
field since it does not have a regular limit as m — > 0, 
where m is the mass. 

As is well known, a theorem of Wald |||!|] plays a cru- 
cial role in this field. The theorem states that if one 
has two different prescriptions for obtaining stress ten- 
sors from metrics and from quantum states, and if these 
prescriptions obey a certain set of physically-motivated 
axioms, then the two prescriptions must agree up to a lo- 
cal conserved tensor ||. Thus, if (T a b) and (T a b) are two 
different such prescriptions for computing the expected 
stress-energy tensor, then the difference 



tab=(T a b) — (T a b) 



(1.2) 



must be a conserved local curvature tensor. 

In the case of a massless scalar field, there is no natural 
mass scale in the theory, and the following well-known ar- 
gument based on dimensional analysis shows that (T ao ) is 
unique up to a two parameter ambiguity. Let us use units 
in which Ti = c = 1, but in which 1. Then, there are 
only two independent conserved local curvature tensors 
with the appropriate dimensions of (mass) 4 , namely 



1 



-9 Sg ab (x) 



-gR(x') 2 (1.3) 



and 



-9 5g ab (x) 



Thus we must have 



t„ 



aH 



(i) 



(1.4) 



(1.5) 



where a and (3 are two unknown dimensionless param- 
eters. Hence, the expected stress tensor (T ao ) is unique 
up to a two parameter ambiguity. 

Consider now a massive scalar field. In this case the 
above argument fails, since there is a preferred mass scale 
present, namely the mass m of the field. Using this mass 
scale one can construct local conserved tensors of the 
form 



1 



-g Sg ab (x) 



d 4 x\ 



-gm 



R(x') 



(1.6) 



which have dimension (mass) 4 and are thus possible can- 
didates for t a b. Here n can be any integer greater than 
2. Similar terms can be constructed from the Ricci and 
Riemann tensors. The conventional view has been to ex- 
clude such terms as being unphysical, since they diverge 
as m — > (see, e.g., p. 90 of Ref. M). Thus, the con- 
ventional view has been that the allowed ambiguity for 
a massive field is no worse than that f or a massless field, 
namely the two parameter ambiguity (1.5). 
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B. Ambiguity in expected stress tensor for a massive 
scalar field 

The first main point of this paper is that the above 
conventional view is unfounded. This can be seen as 
follows. Consider the local conserved tensor 



t a b — 



R(x') 



(1.7) 



where F(x) is any dimensionlcss function of a dimen- 
sionless argument x. In order for t a b to be acceptable on 
physical grounds as a contribution to the expected stress 
tensor, it must satisfy the requirements that 



tab -» 



as 







and 



as R -> 0. 



(1.8) 



(1.9) 



Now, the tensor (1.7) can be written as 

tab = -m 2 R ab F'(R/m 2 ) + m 4 F(R/m 2 )g ab /2 

+F"(R/m 2 )Y ab + F^(R/m 2 )Z ab /m 2 , (1.10) 

where Y a b and Z a b are tensors constructed out of deriva- 
tives of R, of dimension (mass) 4 and (mass) 6 respectively. 
Suppose that we choose the function F to be smooth, 
to satisfy F(0) = 0, and to satisfy F^(x)x j - 2 -> as 
x — > oo for j = 0, 1, 2, 3, where F^ (x) is the jth deriva- 
tive of F. Examples of functions satisfying these require- 
ments are F(x) — x 2 exp(— x 2 ) and F(x) = x 2 /(l + x A ). 
Then, the tensor (1.7) will satisfy the required properties 
(U) and Q ^fj. 

Note that when one expands the function F(x) as a 
power series to obtain 



' — ' L m 



(1.11) 



it can be s een that the tensor (1.7) contains terms of the 
form (1.6), each of which individually has unacceptable 
behavior as m — > 0. However the sum (1.7) of all these 
terms does have acceptable behavior as m — > 0. Note also 



that it is not possible to exclude terms of the form (1.7) 
by postulating, as an additional axiom, that the stress 
tensor be an analytic function of m 2 , since for example 
the choice F(x) = x 2 /(l + x A ) yields a local conserved 
tensor t a b which is an analytic function of m 2 in an open 
neig hborhood of the real axis in the complex m 2 plane 

B 

The ambiguity t a b in the stress-energy tensor allowed 
by the Wald axioms is therefore much worse in the mas- 
sive case than in the massless case. It is an infinite param- 
eter ambiguity — one can specify a free function F(x) — 
rather than a two parameter ambiguity jl2| . Of course, it 
is still possible that the various conventional calculational 



met hod s still agree to within the two parameter ambigu- 
ity ( |l.5| ) |jl3|| . However, there is no guarantee that this 
should be the case. Therefore it would be worthwhile to 
find some additional axiom or physical principle, to aug- 
ment the Wald axioms, that would further pin down the 
stress tensor in the massive case Hi- 



C. Nearly flat spacetimes 

Spacetimes which 
are linear perturbations off Minkowski spacetime form 
a useful testbed in which to probe these issues jD| . The 
second principle purpose of this paper is to explicitly cal- 
culate the renormalizcd stress tensor of a massive scalar 
field in such spacetimes, using the closed-time-path or 
in- in effective action formalism []| [?]]. If our calculation 
is repeated using the point-splitting or deWitt-Schwinger 
methods, then it will be possible to compare the predic- 
tions of the different methods. 

Now the Wald axioms imply that any prescription for 
calculating the stress tensor is determined by specifying 
the expected value of the stress tensor in the incoming 
vacuum state |0,in). The expected value in any other 
state is then uniquely determined Therefore, it suf- 
fices to consider the expected value of the stress tensor in 
the incoming vacuum state. In the massless case, calcula- 
tions of the in- in expected stress tensor have already been 
per formed using several different methods [see Horowitz 
[|l6| and Jordan jnj]]. The results of of these different cal- 
culations agree up to the two parameter ambiguity (ITq), 
as they must according to Wald's theorem. 

The result we obtain fr om t he in-in effective action 
formalism 15H7I is [cf. Eq. below] 



(0, in| f ab (x) |0, in) = aH^ (x) + (3H a l' 0*0 



r( 2 )/ 



1 



256tt 2 



d 4 x' 



H^ix'Wx-x') 
+H a i>(x')T 2 (x-x')\. (1.12) 



Here it is assumed that the metric tensor is of the form 



9ab = T) a b + h a b, 



(1.13) 



where r\ a b is a flat, Minkowski metric, and that the co- 



ordinates x a and 



are Lorentzian coordinates with re- 



spect to rjab- Also a and (3 are arbitrary dimensionless 
constants, the distributions T\(x — x') and T2(x — x') are 



(2) 

Kb 

are linearized versions of the local conserved curva- 



defined by Eqs. ([4. 14[) and (|4.15|) below, and fity and 



ture tensors (1.3) and ( |l.4| ) 



The stress-energy tensor fll.12 ) is causal, as it must 
be, and reduces to the known result of the massless case 
fll6| , [T7"t in the limit m —>■ 0. Furthermore it is not a 



smooth function of 



at m = 0. The calculational 



2 



method we use also automatically yields two undeter- 
mined parameters a and /?, so the result ( 1.12 ) explic- 
itly exhibits the two parameter ambiguity, just as in the 
massless case |T^ |. 



D. Organization of this paper 

Section |H| reviews the in-out and in-in effective action 
formalisms. In Sec. [II we calculate the in-in and in-out 



effective actions of a massive scalar field propagating on 
a spacetime which is a linear perturbation off Minkowski 
spacetime. In Sec. IV we find the expected stress-energy 
in the incoming vacuum state, and then in Sec. |v| we 



discuss its properties. Sec. VI summarizes our results. 

Throughout we use units in which Ti — c — 1, and 
use the metric signature and sign conventions of Misner, 
Thorne and Wheeler pal. Further notational conventions 
are given in AppendixW 



II. THE IN-OUT AND IN-IN PATH INTEGRAL 
FORMALISMS 

In this section we review both the standard in-out 
path i ntegr al formalism of quantum field theory (see, e.g., 
Refs. jlJjU), as applied to curved spacetimes, and the 
modified in-in formalism due to Schwinger 0] . This in-in 
or closed time path method was later adapted to curved 
spacetimes by Jordan |6|Q] , and has been extensively ex- 
plored by H u |2^| . Our presentation of the in-in method 
in Sec. 



[I C below differs from that of Refs. 



in that 

all the fundamental definitions are explicitly coordinate 
independent. 



A. The classical theory 

We consider a massive, minimally coupled scalar field 
for which the Einstein-Klein- Gordon action is 



ah 



S g \9 ab ] + Sm[g ab ,4>], 



where 

and 

S m [g ab , 



S g [g ab ] = 2fx 2 p J J—gR<F 



(2.1) 



(2.2) 



-\ J ^g-{g ab V a <pV b ^ + m 2 <p 2 )d n x. (2.3) 

Here m is the mass of the scalar field, /i 2 — (32ttG)^ 1 is 
the square of the Planck mass, and n is the number of 
spacetime dimensions (we shall be using the dimensional 



regularization scheme below). The corresponding equa- 
tions of motion are Afi 2 G a b = T a b and (□ — m 2 )(f) = 0, 
where the stress-energy tensor is 



-2 S_ 



S m [g ab , 



V Q </)Vb0 - i<? a &V c 0V c - ^g ab m 2 <f> 2 . (2.4) 



B. In-out formalism 

We assume that the metric g a b is asymptotically static 
at early and late times so the incoming and outgoing 
vacuum states |0,in) and |0,out) are well defined. In the 
usual way we define the generating functional 



2 iM/[9 ° b ' J1 = (0,out|0,m)j, 



(2.5) 



where the subscript J on the right hand side indicates 
that a source term 



(J, 4) = / cFxV^gJx) J{x)(t){x) 



(2.6) 

has been added to the action. It can be shown that 
gtW[g a ,J] nag j-jjg integral representation 

e iW[g ab ,J] = J D<f> e l ( s ™[s at >]+W>) 7 (2.7) 

and that time-ordered matrix elements are given by 

{0,ont\T^(x)4>(y)\0,hx)j-- 



xe 



Dcf> <l>(x)4>{y) 

i(S m [g ab ,<t>] + (.J,<t,)) 



(2.8) 



In Eqs. (2.7) and (2.8), the usual boundary conditions on 
the path integral are assumed, namely that <f> is purely 
negative frequency (oc e with lu > 0) at early times and 
positive frequency at late times, or equivalcntly, that the 
mass squared parameter is understood to have a s mal l 
nega tive imaginary part, m 2 — > m 2 — ie. From Eqs. (2.8) 



and (2.4) it follows that 
-2 SW[g ab ,0] 



'-g 



6g 



ah 



= -ie-' lW[9 01 



D(f> 



-g < 



(0,out|f afc |0,in) 
(0,out|0,in) 



(2.9) 



The effective action is defined in the usual way as a 
Legendre transform of the generating functional: 



r m [^]EEiF[ 5 ab ,j]-(j, 



(2.10) 



3 



where 



1 5W[g ab ,J] 



-9 



6,1 



= (O,out|0|O,in)j. (2.11) 



Here and henceforth the subscript m in T m indicates tha t 
the classical action from which T m is computed [Eq. (2.7) 
above] includes only the matter part S m and not the 
gravitational part S g . Combining Eqs. ( |2.9D - (2.11) we 
now obtain 



-2 ST m [g 



ab 



-9 



Sg ab 



<p=cl>[g a 



(0,out|f afc |0,in) 
(0,out|0,in) 



(2.12) 



Here the right hand side is a functional only of the metric 
g ab , while on the left hand side 0[g ab ] is the solution to 
Eq. (|2~Tll) at J = 0. 



Since the action (2.3) is quadratic in <f>, it is straightfor- 
ward to compute the effective action exactly. The result 
is the standard, formal, expression |2(i|] 



ab 



S m [g ab A] 



tr In A 



(2.13) 



where A is the operator given by Atfi = (□ — m + ie)4>. 
The operator A is the natural operator associated with 
the quadratic form S m and with the covariant inner prod- 
uct on functions on spacetime 



(f,h) = / d n x^gjx)f(xyh{x). 



(2.14) 



The reason the inner product (2.14) is the appr opri ate 
inner product is that the measure D(j> in Eq. (2.7) is 
determined by the metric g ab [^l| . 

Finally, we define the quantity T[g ab , 4>] (as opposed 
to r m [g ab ,(j)}) to be the effecti ve a ction obtained when 
one starts fr om t he full action (2J) rather than just the 
matter part (2.3). It is clear that 



r[.9 



a b 



a b 



S g [g 



ab] 



(2.15) 



C. In-in formalism 



We introduce the generating functional 



iW[gi 



,J+,J-] = 



= J Da (0,in\a,T)j_ (a,T\0, 



>J + , 



(2.16) 



which depends on two independent sources J+(x) and 
J-(x), as well as two metrics g^ and g°^. Equation 



(2.16) includes an integration over a complete set of 



states | a, T) on the hypersurface x° = T at some fu- 
ture time T. We assume that gf = g°L b for x° > T. Each 



of the matrix elements in Eq. (|2.16 ) can be expressed as 
path integrals in the usual way: 



(o,r|0,in) J± = / Dfee* 3 ^ ,<P±]+(J±,<I>±)±) ^ ( 2 .17) 



where the measure D(j) + and inner product (. . .) + are 
determined by the metric <?" b , the measure D(j>_ and the 
inner produ ct (. . .)_ by the metric g°^ ' . Combining Eqs. 
( pT6| ) and ( p7l7| ) yields 



iW[g a , 



,J+,J-]- 



= / Da 



D(j)_e 



-i(S' m lg! b ,<P-] + (J-.,4>-)-) 



D(j)+e 



with the boundary condition that 



i(Sm[af ,<!>+]+ (J+,4>+)+) 



(2.18) 



a on the 



hypersurface given by x° — T. Another boundary condi- 
tion, needed to assure convergence of the path integrals, 
is that (/)+ be purely negative frequency and 4>- be purely 
negative frequency at early times, or equi valc ntly that m 2 
be interpreted as m 2 — ie in the action (2.3). From now 
on we assume the second of these. We can rewrite the 



generating functional ( 2.1£ ) as 

p W[gf ,g a _ b ,J+J-] 

D<t>+D<l>_expi{s m [gt, 4>+] + {J+;4>+) + 
-S* m [gl b ,cl>_]-{J_,^)_}, 

(2.19) 

where the integral J Da is now included in the integra- 
tion over cj) + and 4>- , and the boundary condition is that 
4> + = tfi_ on the hypersurface given by x° = T. Bel ow 
we will be taking the limit T — > oo. From Eq. ( 2.19 ) it 
follows that 



-2 6W[gf, g a _ b ,0,0] 



Sgf 



ab 



gl-g_=g a ° 

= (0,m|f o6 |0,in). (2.20) 



The effective action is defined to be a Legendre trans- 
form of the generating functional as before: 



r 



m [9+ ,9- ,«?+)« 



(2.21) 



; W[g+, g° b , J+, J_] — (J + , 4>+)+ + (J-,4>-)-, 



where 



i±[gt,g ab ,J+,J-}=± 



1 5W[gt,g ab ,J + ,J-] 



5J± 



(2.22) 



and where we use the shorthand notation (2.6). Note 
t hat when g±_= g ^ and J + = J_, we find from Eqs. 
( pTl9| ) and ( p72l| ) that 

4>±[g ab ,g ab ,J,J} = (0,hM0,m).,. (2.23) 



4 



From Eqs. fl2.20| ) and fl2.2l| ) we obtain the expected 
stress-energy tensor in the in-coming vacuum state: 



b in— i 



(0,in|T o6 |0,in) 
-2 6T r , 



(2.24) 



Here 4>[g ab ] is given by Eq. (2.22) at J+ = J_ = and 



„afc „ a 6 ^ab. 



(2.25) 



From Eq. (2.23) it can be seen that 4>[g ab ) is the expected 
value of the field in the incoming vacuum state. Equation 
( p. 24 ) is the formula we will use below to compute the 
stress tensor. 

We introduce the shorthand notations 6* = 



J, = 



-J+ 
-J_ 



(2.26) 



and 



S^f^+j-S^,^}, (2.27) 



where the index s take s the values + and — . The gener- 
ating functional (2.19) can be rewritten using these no- 
tations as 



JW[gf ,J,] _ 



Dcj) r e 



i(S m [gf,4 s ]+(Jt,4, t )t) 



(2.28) 



where a sum over the repeated index t is understood. 

Nex t, we derive the analog of the formal expression 
(2.13) for the in-in effective action. We define the oper- 
ator A on pairs of functions <j) s = (<f) + ,4>-) by 



(A(f>) s (x) = J d n yy/-g t (y)A st {x,y) </> t (v), 
where 

5 2 S m [gf^ r ] 



1 



y/g a (x)gt{y) H s (x)Ht(y) 



(2.30) 



which can be written using Eqs. ( |2.3| ) and (2.27) as 

A st (x,y) = S'\x-y) 



D + ,x — m +ie 









(2.31) 



Here □+ x denotes the wave operator associated with the 
metric g^ acting on the coordinates x — x a , and simi- 
larly for a_ tX , Using Eqs. ( |2.21[ ), (|2.28|) a nd ( |2.30| ), we 
then find the following analog of Eq. (2.13) 



T m [gf, 4> r ] = S m [gf, 4> r ] + -Tr In (A st ) , (2.32) 



where Tr denotes the appropriate trace over both space- 
time variables and over the indices s, t. The operator A is 
t he n atural operator associated with the quadratic form 
( |2.27|) and with the covariant inner product on pairs of 
functions (/+, /_) given by 

((f+J-),(K,h-))= Yl l d n x^gAx)f s {x)*h s {x). 



s=+, 



(2.33) 



The reason the inner product (2.33) is the appropriate 
inner pro duct is that the measures D<f> + and in 
Eq. ( [2.19] ) are determined by the metrics and g ab 
respectively. 

In our perturbative computations below, we shall de- 
rive an expression for the effective action in terms of a 
series of products of operators. For that purpose, it will 
be convenient to use the Hilbert space structure associ- 
ated with the coordinate dependent inner product 

((/+,/-), (/H,M) C = E [ d n xf s ( X yh s (x), 

s=+ - J 

(2.34) 

instead of that associated with the covariant inner prod- 
uct (2.33). We will a lway s choose the coordinate system 
x a appearing in Eq. (2.34) to be a Lorentzian coordinate 
system associated with the flat metric rj a b- 

Finally, we define the quantity r[g" b , g ab , (f> + , </>_] (as 
opposed to T m ) to be the effective action obtained when 
one starts fr om t he full action (2.1) rather than just the 
matter part (|2.3|). It is clear that 



r[gf , g a _ b 7 4> + , 0_] = r m [gf , g a _ b , 0+, £_] 

+S g [gt]-S g [g ab }. (2.35) 



(2.29) The semiclassical equations of motion are given by 



ST 



Sg 



a b 



9T 



ST 

Hi 



0. (2.36) 



91" 



The second of these equations is a utom atically solved 
when we choose cj> = (j)[g a % cf. Eq. (gjgj) above, corre- 
sponding to the incoming vacuum state. Thus the equa- 
tion of motion reduces to 



ST 



Sgf 



(2.37) 



4>±=4>[g ab } , g±=g a 



III. THE EFFECTIVE ACTION FOR NEARLY 
FLAT SPACETIMES 

In this s ectio n we specialize to almost flat spacetimes of 
the form (1.13), and calculate the in-out effective action 



5 



( [2.15 ) and the in-in effective action (2.35) as series expan- 
sions in powers of the metric perturbation h a b- We use 
the methods of Hartle and Horowitz pj| and of Jordan 
P,^7| , who performed similar calculations in the massless 
case. 

We start by further simplifying Eq. ( 2.32| ). Note that 



Sg 



- b Trln[A st ] = Tr (A 



Sgf 



tr A 



_! SA 



++ 



Sg 



ab 



(3.1) 
(3.2) 



since in the sum over r in Eq. ( |3.1| ), only A ++ depends 
on gf. Here tr denotes a trace over the spacetime vari- 
ables x, y only, n ot inc ludin g a s um over the s, t indices. 
Combining Eqs. and @ we now obtain 



-t,\nA ++ +F 1 [g ab ], (3.3) 



where F\ [<?" b ] is some functional of g_. This term F\ [g ^] 
will n ot contribute to the functional derivative in Eq. 
( [2.24 ) and hence will not contribute to the in-in expected 
stress tensor. 

Next, we define the coordinate dependent operator 
A = A rs acting on pairs of functions (<^+ , <p- ) by 



{Acj>) r {x) = ^g7^){Act>) r {x) 



(3.4) 



Then from Eqs. J2.29| ) and ( [2~3l] ) the kernel A rs (x,y) of 
A wit h respect to the coordinate dependent inner product 
( [2.34 ) is given by 



A rs (x,y) = y / -g r (x)y / -g s {y)A rs (x,y). 



(3.5) 



[Here A rs (x,y) is the kernel of t he operator A w ith r c- 
spect to the inner product given by Eq. (|2.3l|) 1 

Combining Eqs. ( fOl] ) and (p.3[) — (Q now yields 



T m [gf, 4> r ] = S m {g? b , r ] + Fxig^ 



+- tr In A, 



F 2 [det(<7 + )], 



(3.6) 



where -F2[det(<?-|_)] is some functional of the determinant 
det[gf(x)}. When we take the variational derivative of 
r ?n in Eq. ( ]2.24[ ) to calculate the stress tensor, the term 
F2 will contribute a term proportional to the metric g a b 
and hence will contribute only to the renormalization of 
the cosmological constant. 

We define the propagator G s t to be the inverse of the 
operator A st . Its kernel G s t(x,y) with respect to the 
inner product ( |2.33 ) is given by 



(3.7) 

Note that the operation of taking the inverse is unique by 
virtue of the boundary conditions imposed on the path 



integrals and the fact that the mass squared parameter 
is assumed to have a small negative imaginary part (see 
Appendix ^). We also define a coordinate dependent 
operator G s t to be the inverse of the operator A st ; it s 
kernel G s t(x,y) with respect to the inner product ( 2.34 ) 
is given by 



d n x' A rs (x, x')G st (x', y) = -5 n (x - y) S rt . (3. 



From Eqs. J3J), (fcj]) and (^8|) it follows that 
G rs (x,y) = G rs (x,y), that is, the above two kernels co- 
incide. 



A. Perturbation expansion for the in-in effective 
action 

We now expand the operator A st as 

A rs ( X) x') = A° s (x,x') + V rs (x,x'), (3.9) 

where 

V rs (x,x') = V}p(x,x') + V r f(x,x') + ... (3.10) 

Here A® s is the Minkowski spacetime operator, and the 
terms and are the pieces of A rs that are lin- 
ear and quadrati c in t he me tric p erturbation h a b, respec- 
tively [see Eq^l^l) and ([3.28 ) below]. Note that from 
Eqs. ( 2.31 ), (3.5) and (3J;) it follows that the operator 
V rs is diagonal in the indices r and s and is of the form 



ab „ab~\ 



v rs [g+,g 



V ++ [gf,g a _ b ] V + -[gf,g ab ] 
V- + [gf,g ab ] V—[g ab ,g ab ] 



V[gf} 






-V[g ab }* 



(3.11) 



for some functional V — V[g ab ]. 

We similarly expand the propagator G st as 



Gst — G° st + G 



(i) 



G 



(2) 



(3.12) 



The Minkowski spacetime pro pa gato r G" t can be ob- 
tained by combining Eqs. ( 2.31 ), ( ft.5| ) and (3.8) and us- 
ing g a b — rjab- In Appendix ^| we show that this yields 



G° st (x,y) 



where 



G{x) 



G° ++ (x,y) G° + _(x,y) 
G°_ + (x,y) G°__(x, y) 



G{x~y) A+(x-y) 
-A+{x-y)* -G(x-y)* 



d n p e lpx 
(27r)™ p 2 + m 2 — ie 



(3.13) 



(3.14) 



G 



is the free Feynman propagator and 
d n p 



A + (x) = 2ni 



lpx 5(p 2 + m 2 )0(~p°) (3.15) 



effective action (3.17) and using Eqs. ( fUll ) and ( |3| ), 
we see that the only difference is the term involving 
G°_ + (x, y), since G° + = G° and V+ + = V. Hence we 
can write 



is the positive Wightman function. 

Next, by combining the e xpa nsions (3.E) and ( 3.12 ) 
together with the definition (3.S), we find the following 
expression for the logarithmic term appearing in the ef- 
fective action (3.3) (see Appendix ^): 



ln[G ++ ] = ln[G° ] + V^G 



+0(h 3 ). (3.16) 



The products on the right hand side of Eq. ( |3.16 ) are 
operator products, where t he ke rnels are understood to 
refer to the inner product ( |2.34| ), so that, for example, 

\y$GP ++ Ux 3 v) = Jp/vl¥&y')G° ++ (y',y). Combin- 
ing Eqs. ( 3.16 ) and ( |3.3|) we finally obtain the perturba- 
tion expansion of the in-in effective action 



-tr \nG°, , - -tr 
2 ++ 2 



1/(1)^0 



1 



2 v ++^++ ' 



+V«G° 
+F 1 [ 5 ° b ]+F 2 [det( 3+ )] + 0(/i 3 ). 



(3.17) 



r m [«4 6 ,5- b ,0+,<Min-m 



where 



r m K h ,0+]in-out+F 3 [.g^] 

+l%f,^ fc ,0 + ,0_], (3.19) 



U[gf 1 g ab J + ,^] = --tr 



^G^_^G°. 



(3.20) 



and Fs\g_ ] is a term which does not depend on g^r or 

C. Explicit calculations 

We write the spacetime metric as 

9ab = Vab + Kb, (3-21) 

where rj ab is a flat Minkowski metric. From now on in- 
dices are raised and lowered with ij ab , and derivatives 
denoted by a comma are coordinate derivatives in a 
Lorentzian coordinate system associated with the met- 
ric r\ ab . Expanding the action (2.1) to second order in 
h ab yields (see Appendix [a|) 

S[g ab ^]^S g [h ab ]+S m [h ab ^]+0{h 3 ), (3.22) 



where the quadratic actions S g and S m are given by 



B. Perturbation expansion for the in-out effective 
action 

Consider now the corresponding calculation in the in- 
out formalism. If we compare the gener ating functionals 
(2.7) of the in-out formalism and ( 2.19] ) of the in-in for- 
malism, we see that the terms in <j) and in <f> + in these 
equa tions coincide. Hence, from the definitions (2.1C ) 
and ( 2.21 ) of the effective actions in terms of the generat- 
ing functionals, it follows that the in-out propagator and 
effective action can be obtained from the corresponding 
in-in quantities simply by replacing (f> + with </>, J + with 
J, and dropping all terms generated by and J_ . The 
resulting in-out action is 

T m [g ab ^} = S m [g ab ,4>] - -tr lnG° - l -tr[v^G° 

+ y(2) G o + l y (i) G o y (i) G oj + F2 [ det(5)] 
+0(h 3 ), (3.18) 
where G°(x, y) = G(x — y) is the usual free Feynman 



propagator (3.14), and V' 1 ) and are the pieces of 
the operator V defined in Eq. ( |3.11) tha t are linear and 
quadratic in h ab . Comparing Eq. ( |3.18| ) with the in-in 



S g [hab]= 2,4 / (h ab ' ab - h a a )d n 



(3.23) 



and 

S m \h a bi ' 



(V ab ^,b + m^ 2 )d n x 



~fo?°fy«0,6 - h ah <t>, a ct>, b + hm 2 ^ 

(I A 3 _ h ab h ab ) V ab <P,a(t>,b 
a 4 



+ (h 2 - \h ab h ab )m 2 <f - \hh ah (j>, a ct>, b 

o 4 Z 



+h ac h b c cp^ fi 



d n x. 



(3.24) 



Here h is the trace h = rj ab h g b . 

Next, we find from Eqs. ( |2.31 ) and ( p\5| ) the formula 
for the operator A ++ 

A ++ (x, y) = 5 n (x - y) ( V ab d a d b - m 2 + V), (3.25) 



where [cf. Eq. ( 3.11 ) above] 

V = V^ +V^+Q(h 3 ) 



(3.26) 
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with 



yW = -d a hd a - d a h ab d b - -m 2 h 



-d a h ab d b 



i 2 h 



(3.27) 



and 



V® = d a h ac h b c d b - \d a h ed h ed rf b d b - — 
4 n — 



-d n hh ah d h 



1 



4(n - 2) 

2 



+ - 



A// 



h cd h. 



d a h 2 r) a % 
1 



(3.28) 



Here for simplicity we have written h a b+ simply as h ab . 
Also we are using a notational convention where, for ex- 
ample, (d a hd a )cp = d a (hd a tp). We have also introduced 
the quantity 



1 



-hl]ab 



(3.29) 



and its trace h = n ab h ab . Note that h ab is the trace- 
reversal of the metric perturbation h ab when the dimen- 
sion n of spac etime is 4, b ut not ot herw ise. 

Using Eqs. (f]||, ( |3~ll"l) , ( |3~27| ), and we 

now find that the in-in effective action is 



T[h a b+,h a b-,i>+, <t>-] = S[h ab+ ,h ab -,(f> + ,4>-] 
+U[h ab+ ,h a b-} + Ki[h ab+ ] + K 2 [h ab+ } 
+L[h ab+ ] + F^g 11 }, det(<7+)] + 0(h 3 ), 

where 



(3.30) 



S[h ab+ ,h ab _,4> + ,(j)_]= S g [h ab+ \ + S m [h ab+ ,<f> + ] 

-S g [h ab -\- S* m [h ab -,(f>-]. (3.31) 

In Eq. ( |3.30| ), we have absorbed the constant tr In G° ++ 
and the terms F% and F 2 into the functional 
Fi[g ab , det(g+)]. As before this term is a quantity which 
depends on the metric g°^ only through its determinant, 
and which thus affects the in-in equation of motion only 
via a renormalization of the cosmological constant. The 
term s U\h ab+ , h ab J\, Ki[h ab+ ], K 2 [h ab+ ] and L[h ab+ ] in 
Eq. (3.30) are given by 



U[h ab+ ,h ab _] =- % -j d n x j d n x'V^[h ab+ {x)\ 

xG° + _(x - x') {-V^>[h ab -{x')]) G°_ + (x' - x), 

(3.32) 



K 1 [h ab ]=-~ j d n x / d n x'S n (x-x')V w (x)G(x-x') 



d n x / d n x'8 n (x-x') 



d a h ab d b 



n-2 



i 2 h 



G(x-x'), 



(3.33) 



K 2 [h ab ] = -- / d n x / d n x'S n (x - x')V {2) (x)G(x - x') 



d n x / d n x'S n (x - x') 



d a h ac h b c d, 



X -d a h cd h cd if b d b - -^—d a hh ab d b 
4 n — 2 



1 



d a hhrj ab d b 



4(n - 2) 



■{h cd 'h. 



cd 



n-2 



hh) 



G(x - x'), (3.34) 



and 



L[h ab ] = -- I d n x I d n x 



x V<M{x)G(x - x')V {l \x')G{x' - x) 
= Li [h a b] + L 2 [h ab ] + L 3 [hab] + Li [h ab ] . (3.35) 



We have also defined 

i 



Li[h a 



d n x / d n x' 



xh(x) ab G(x - x').b c> h{x') c ' d 'G{x' - x), d -a , 

(3.36) 



L 2 [h ab ] =— I <V.r I r/'V 



x -^—h(x)G(x - x'). c ,h(x') c ' d 'G(x' - x), d , , 
n — 2 

(3.37) 



L,\h 



3[ n ab\ 



d n x / d n x' 



x h(x) ab G(x - x'), b -h(x')G(x' - x), a 



n-2 



(3.38) 



and 

Li[h ab ] 



--J d n x J (Fx' 

2 2 
777 _ 777 _ 

x -h(x)G(x - x') -h(x')G(x' - x) 

n — 2 n — 2 



(3.39) 



In Eqs. (3.32) and (3.35), the differential operators in 
each factor of act only on the propagator immedi- 
ately to the right of such factors. 

It is straightforward to obtain the in-ou t effective 
action from th e in- in effective action (|3.30[) usin g the 
method of Sec. |IIIB| above. Using Eqs. |TTg) and ( |3.30| ) 
we obtain 
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r[h ab 



S[h ab , 0] + K^hab] + K 2 [h ab ] + L[h ab ] 
+F 2 [det{g)] + 0(h 3 ). (3.40) 



As already mentione d [cf. Eq. ( |3.19| ) above] the term 
U[hab+,h a b-] m Eq. ( 3.32 ) contains the differences be- 
tween the in-in and in-out formalisms, and has to be 
added to the in-out effective action T[h a b, (f>] to obtain 
the in-in effective action T[h a b+, <t >+ » fap ,fe- 1 4>-]- 

Next, we insert the expressio n (3.14 ) for t he Feynman 
propagator G(x — x') into Eqs. (3.33) — ( 3. 39] ) to evaluate 
the quantities U, K\, K 2 and L. To simplify the calcu- 
lations, we work in the Lorentz gauge where h ab b = 0, 
and we regularize the results using dimensional regular- 
ization. The results are written in terms of curvature 
invariants using Appendix |b], and are listed in Appendix 



As an example, we show how to compute the term 
(HI): 



Lz[h a b] = 



d n x / d n x'h a \x)- 



4 

d n p 

(2tt)™ J (2tt)™ (p 2 + m 2 - ie)(q 2 + m 2 - ie) 



Pbq a e 



-h(x') 

i{p-q)(x-x') 



- \ d n x / d n x'h{x) al 



d n k 



m ~ 



:) ik(x—x') 



n-2 
d n q 



h(x') 



(27t)" J (27r) r 

q a (h + qb) 



((fc + q) 2 + m 2 — ie){q 2 + m 2 — ie) 



(3.41) 



We can drop the kb in the term {kb + qb), since we are 
working in the Lorentz gauge. This yields 



L 3 [hab] = — A \ d n x / d n x'h{x) a 



-Ux') 



d n k 
(2ir) n 



(3.42) 



where 



q a qb 



((k + q) 2 + m 2 — ie)(q 2 + ' 



ie) 



(3.43) 



In order to perform this integral, we analytically con- 
tinue to Euclidean signature. Defining /^(fc ,^ 3 ) = 
I a b(—ik°, fc J ) and using the transformations 



q° -> -iq° and fc° -ik° 
in Eq. ( |3.43| ), we obtain 

id n q 



(3.44) 



Iabik) 



dx 



q a qb 



[q 2 +m 2 + 2qk(l - x) + k 2 (l - x)} 2 



(3.45) 



In Eq. (3.45) we have also introduced the Feynman pa- 
rameter x, and dropped the ie. The integral over q can 
now be evaluated (see, e.g., Ramond (24|]). The result is 



Iabik) 



dX (4 7 r)«/2r(2) 



+ 



(t) 5 °+ s °6 ab T(l-n/2)/2 
[m 2 + k 2 (l -x)- k 2 il - cc) 2 ] 1 -"/ 2 



(3.46) 



Next, the term k a kbf(k,x) can be dropped, since it 
will not contribute to Lz[h a b\ as we are working in the 
Lorentz gauge. Also, when we analytically continue 
back to Lorentzian signature using k° — ► ik , we find 



(0 



hh{k) 



Tjab- 



Thus we obtain 



Vab 



(4tt) 2 2 - n 



— 7 + In 47T 



( / dx [in + k (1 - x)x] 
'o 

X |l - (2 - n/2) ln[m 2 + fc 2 (l - x)x]\ 
+0[(4-n)], (3.47) 

where 7 is Euler's constant. The integral over x can now 
be performed with the result 



hb(k) 



Vab 



(4tt) 2 2 - n 



4 - n 

V 6/ 3 3 V 



7 



m 

47T 
2 1 /i™2\ 



4m 2 ) 



fc 2 + 4to 2 
k 2 

0[(4-n)]. 



arctanh 



fc 2 



fc 2 + 4m 2 



(3.48) 

Below we will write the function arctanh(if ) = ln[(l + 
K)/(l — K)]/2 appearing in Eq. ( 3.48| ) in terms of a log- 
arithm. When the ie from the mass-term is included this 
will lead to logarithms of the form 



ln(K - ie)= ln \K\ - iirQ(-K) 



(3.49) 



where Q(K) is the Heavyside step function. Henceforth 
whe n we write ln we shall mean the logarithm defined in 
Eq. ( 3.49| ), which has a branchcut along the negative real 
axis. 



When the remaining terms in Eq. (3.30) are evaluated 
in the same fashion, and written in terms of curvature 
invariants (see Appendix |b|) , we find the results listed in 
Appendix pi. The effective action (3.4C) then becomes 

r[h ab ,4>] = s[h ab ,4>} + w, (3.50) 

with 
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w = 



1 



5127T 2 

. Coo 



d 4 x 



8m 4 A,o V^g - 8— Boo V^gR 



+W nl + F 2 [det{g)]. 



(3.51) 



Here 



1 



d 4 x I d? x > I ALrW*-*') 



512tt 2 J J J (2?r) 4 

i?(x)i?(x')Qi(fc) + R ab (x)R ab (x')Q 2 (k) 

(3.52) 

is the non-local part of the effective action. Furthermore 
we have defined the functions 



L (fc)=4 



4 to 37 to 2 1 to 2 — ie 
In ■ 



15 fc 4 45 fc 2 30 ^ 2 
(fc 2 + 4w 2 ) 2 2m 2 (fc 2 +4m 2 ) 2m 4 " 



30fc 4 



3k 4 



4(m 2 -ie), /v 1 
xi 1 + — — — -m 



k 4 I 

4(m 2 — ie) , -i s 
P h 1 



fc 2 



4(m 2 — ie) 



k 2 



1, 



(3.53) 



and 



Q 2 (fc) 



16 to 4 28 m 2 1 m 2 — ie 
15 fc 4 " + 45 "fc 2 " ~ 30 n /I 2 



30fc 4 



(fc 2 + 4m 2 ) 2 i 4(m 2 - ie) 



x In 



4(m — zg) 



1 



4(m 2 — ie) 

— p — 



- U J 



(3.54) 



The logarithm used here is the o ne de fined in Eq. ( p.49[ ) 
The constants appearing in Eq. (3.51) are 



Aoc = h ln47T - 7 - ln/i 2 + ~ - In 

4 - n 2 /i- 



Note that both Qi(fc) and Q2(k) are finite at fc 2 = and 
that they reduce to 



Q 2 (fc) = 2Qi(fe) = -Aki 



15 

4 r 

15 



In 



-i7r6(-fc 2 ) (3.59) 



if to — 0. Not e als o that the constant /z 2 appearing in 
Eqs. ( 3.53 ) — ( 3.58 ) drops out w hen these equations are 
inserted in Eqs. ( |3.51 ) and (3.52) in the to ^ case. The 
constant \j? has dimension (mass) 2 and has been inserted 
t o yie ld the corr ect dimensions in the logarithms of Eqs. 
( |3~53| ) and (|3~54l ). 



Now the functional i^det^)] in Eq. ( 3.51 ) must be a 
coordinate invariant, since T m and the rest of the terms 
in that equation are. It follows that 



F 2 [det( 5 )] oc / d n xy/^g{x) 



(3.60) 



D. Renormalization of the in-out effective action 



We now rewrite the classical action ( |2.l| ) in terms of 



some bare coupling constants /i 2 A&, aj, and j3\ 



+ / d n xy^g 



2 M 2 b ( J R-2A b ) 



(3.61) 



From Eqs. (|3.50|) , ( |3.5l| ) and ( |3.6l|) we then find 
r[/i„ 6 , 0] = ~ J d 4 x^—g [V a 4>V a 4> + m 2 4> 2 ] 



2/i 2 (i?-2A) 



c? x\/—g 
-(3R a bR ah — ~^°<-R 2 



+VF„i 



(3.62) 



4- n 



ln47r — 7 — ln/i 2 + 1 — In — , (3.56) Here /i 2 , A, a and /3 are the renormalized values of the 

A* parameters, given by 



and 



Coo = -r— — + In 4tt - 7 - In A i 2 + -| (3.57) 
4-n 15 



Boo = t^— +ln47r-7-ln/i 2 + 4- (3-58) 
4 — n 15 



384^ m2B °°' 



; m 4 A c 



256?r 2 



(3.63) 



(3.64) 
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and 



a = a b 



19207r 2 ^°° ! 



,C 



(3.65) 



(3.66) 



960tt 2 

In the usual way, the renormalized values of the parame- 
ters are finite whe n we choose the bare parameters suit- 
ably. In Eq. ( 3.64 ), the quantity A is the (uncalculated) 
contribution to the renor maliza tion of the cosmological 
constant due to the term ( 3.6C| ). Note also that the pa- 
rameters a and (3 which appear in the local part of the 
renormalized effective action and the parameter /i which 
appears in the n on-lo cal part W n \ are n ot all indepen- 
dent: from Eqs. ( 3.52| ) - ( 3.54 ) and ( 3.62| ) it can be seen 
that a change in [i can be compensated for by changes in 
a and (3. 

Finally, the renor maliz e d in-in effective actio n is given 
by combining Eqs. ( |3.19D , fl3.62| ), and Eq. (|e|) from Ap- 
pendix ||] below. 



In writing this tensor we have also introduced the lin- 
earized versions H^(x) and (x) [see Eqs. ( B14 ) and 
(B15) below] of the conserved local curvature tensors 



1 



-9 5g ab (x) 



2g ab UR - 2V a V b R - -g ab R 2 + 2RR ab (4.6) 



and 



-9 5g ab {x) 



(Tx'^R cd {x')R cd {x') 



1 



-g ab UR + DR ab - 2V c V a i? c 



-2R c a R cb — —gabRcd I ' 



cd 



(4.7) 



As is well known, Eq. (4.4) is not a physically realistic 
equation for semiclassical gravity since the right hand 
side is complex and not real PJl7j]. 



IV. THE STRESS ENERGY TENSOR 

A. Equations of motion in the in-out formalism 

The semiclassical equations of motion are obtained 
from 

-2 S 



and 



-g $g ab 
5 



0. 



T[h cd , <P] = 0. 



(4.2) 



Equation (4.2) is automatically solved when we choose 

t = t[g a % (4.3) 

the functional g iven by Eq. J2.ll ) above a t J = 0. When 
we insert Eqs. fl3.62| ) and (^.3|) into Eq. ( p~l| ), the func- 
tional derivative of the first term on the right hand side 
of Eq. (3.62) ca n be dropped since it is of order 0(h 2 ) 
[as 4> from Eq. fl2.1l| ) for J = is of order 0(h)]. The 
resulting equation of motion is 

G a b{x) + Ag ab (x) = -^T ab (x)- m _ out , (4.4) 

where we have defined the in-out expected stress-energy 
tensor 



T a b(x) 



(0,out|f ab |0,i 
(0,out|0,in) 

1 f d 4 x' 



d 4 k 



ik(x—x') 



256tt 2 ./""./ (2tt^ 4 



(4.5) 



B. The in-in expected stress-energy tensor 



By combining Eqs. ( p^ ), ( p7| ), ( p0| ), ( pC| ), Q 
and (4.4) we obtain the equations of motion in the in-in 
formalism 



(4.1) G ab {x) + Ag ab (x) = [T Qb (x) in _ out + T' ah {x)] , (4.8) 



4/y 



where T ab (x)i n - ou t is given by Eq. (4J3), and whe re the 
additional term T' ab (x) due to the term U in Eq. (3.3C) 
is given by 



T' ab (x) 



-2 5U{h ab+ ,h ab _} 



V=9 &9+ 

-1 [ fi*-' f rf4fc ik(x-x') 
256^ 2 J J (2tt) 4 



X 



(4.9) 



Here we used Eq. (ES) from Appendix [e] below, and have 
defined 



Q[(k) = -4 



and 
Q' 2 (k) 



(k 2 + Am 2 ) 2 2m 2 (k 2 + 4m 2 ) 2m 4 



30k 4 



3fc 4 



it 4 



4m 2 



x\ l + — r 2irie{-k 2 -Am 2 )e{-k°) (4.10) 



(fc 2 +4m 2 ) 2 
30/c 4 



l + ^-2mQ(-k 2 -4m 2 )e{-k°). (4.11) 



n 



The in-in expected stress-energy tensor 
_ (0,m|T ab |0,m) 



is therefore given by 

TofoyXjixx—i 



(0,in|0, 



(4.12) 



Tab{x)'ui — out + T'ab(x) 

-l r , 4 , /■ d 4 fc 

ax 



,ik{x — x ) 



256tt 2 7 7 (2tt) 4 



Here we have defined 



h-(2), 



Ti(fc) =Qi(fc)+Qi(A) 

_ f 4 m 4 37 m 2 



15 A: 4 45 k 2 30 



1 TO 

In — 



(k 2 + im 2 ) 2 2m 2 (k 2 + 4m 2 ) + 2m 4 



30k 4 



3fc 4 



fc 4 



/ 4(to 2 — ie) 



x In 




1 4( m 2-» e Sgnfc") -. v 

1 H p + 1 



4(m 2 -t6Sgnfc°) 



- 1, 



(4.14) 



and 

f 2 (k)= Q 2 (k) + Q' 2 (k) 

16 to 4 28 to 2 
15 fc 4 " + 45 'k 2 ~ 30 

(k 2 + Am 2 ) 2 



In- 



30fc 4 



4(m 2 — ie) 



1 



A- 2 



x In 



v / 1+ 4(m'-«S^fc°T +r 
/-^ ^ 4(m 2 — ieSgn/c ) ^ 



(4.15) 



It is easy to see that Ti(fc) and T 2 (k) are finite at k 2 = 
(for to 7^ 0), and that they are sufficiently regular that 
their Fourier transforms T±(x) and T 2 {x) exist as distri- 
butions. 



V. PROPERTIES OF THE IN-IN EXPECTED 
STRESS-ENERGY TENSOR 



The stress-energy tensor given in Eq. ( 4.13| ) is deter- 
mined by t he Gr een functions T±(k) and T 2 (k) in Eqs. 
(4.14) and (4.15). In this section we show that in the 



limit to — > 0, these Green functions reduce to the previ- 
ously obtained Green functions for a massless field. We 
also show that they are causal, i.e., that their Fourier 
transforms T\(x) and T 2 (x) have support only inside the 
past light cone. These properties serve as a check of our 
calculation. 



A. The massless limit 



The Green functions Ti(fc) and T 2 (k) in Eqs. ( |4.14| ) 
and (4.15) reduce to 



(4.13) f 2 (k) = 2f x (fc) 



_4_ (k 2 - ieSgn(fc°) 
15 

4 
15 



In - i 7re(-fc 2 )Sgn(fc°) 

(5.1) 



if to = 0. This Green function together with Eq. ( 4.13| ) 
yields exactly the same the stress-energy tensor as found 
by Horowitz Jll| and Jordan [fl7jl . 

Note that the Green functions are not smooth in to 2 
near to = 0. For fc° > we find 



dT 2 



dm 2 



m 2 =0 



28 
45k 2 



1 



T5fc 2 ,/1-^ 



1 1 / I -|T 
15k 2 



'1 - 



x In 



'1 - 



1, 



(5.2) 



which diverges in the limit e — > 0. Hence the first deriva- 
tive of the stress tensor with respect to to 2 does not exist 
at to = 0. 



B. Causality 

It is difficult to find the Fourier transforms Tj(x), j — 
1, 2, of the Green functions (4.14) and (4.15). However, it 
is not necessary to explicitly perform these Fourier trans- 
forms in order to demonstrate causality. By Lorentz in- 
variance it is sufficient to show that 



and 



Tj(t, 0,0,0) = for t < 



Tj(0,r,0,0) = for r ^ 0, 



(5.3) 



(5.4) 



for j = 1,2. In other words, the Green functions Tj{x) 
must be zero inside the past light cone and outside the 
light cone. To check the condition ( |5.3|) we write 



^(t, 0,0,0) = 



(2^) 4 



d A k / dk»e~ lk %{k). (5.5) 
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Now from Eqs. (4.14) and (4.15) we see that the loga- 
rithmic terms in both T\(k) and ^(/c) have branchcuts 
in the lower complex k° plane, but no poles elsewhere. 
It is therefore possible to deform the contour of the k° 
integration into the usual semi-circle with infinite radius 
in the upper complex A; plane. Since t < the integral 
vanishes. This immediately shows that Tj(t, 0,0, 0) = 
for t < 0. A similar argument can be used to show that 
2K0,r,0,0) = 0. 



VI. CONCLUSIONS 

We have shown that the Wald axioms determine the 
stress-energy tensor (up to two parameters) only in the 
case of a massless field. In the case of a massive scalar 
field, the Wald axioms allow for a much larger ambiguity. 
We have calculated the expectation value of the stress- 
energy tensor in the incoming vacuum state for a massive 
scalar field on any spacetime which is a linear pertur- 
bation off Minkowski spacetime, generalizing an earlier 
formula of Horowitz Jl(| and Jordan j^| in the massless 
case. In our calculation we used the in-in effective action 
formalism [||-f7) . As expected, the resulting stress-energy 
tensor is causal and reduces to the known result in the 
massless case in the limit m — * 0. As in the massless case, 
we find a two parameter ambiguity in the stress-energy 
tensor, even though this is not guaranteed by the Wald 
axioms. 

After this paper was submitted for publication, we be- 
came aware of related work by Dalvit and Mazzitelli [p5| . 
Dalvit and Mazzitelli calculate an in-in effective action 
using the method of expanding in powers of curvature 
pion eered by Vilkovisky and collaborators p^| . Our re- 
sult (1.12) can be derived from Eqs. (38) and (41) of Ref. 
ps| by evaluating the integrals over t, by specializing to 
linear perturbations about flat spacetime, and by repre- 
sentating the d'Alembcrtian operators in Fourier space. 
Our derivation of the result ( 1.12Q has the advantage that 
it is more direct and simple than that of Ref. |25| , since 
it does not rely on the Vilkovisky expansion formalism. 

We conclude by listing some open questions. First, are 
there additional axioms which would reduce the ambi- 
guity in the stress-en ergy tensor? Second, will the same 
stress-energy tensor ( 1.12 ) be predicted (up to the two 
parameter ambiguity) by the point splitting method || 
or by the deWitt-Schwinger method Q? There is no a- 
priori guarantee that this will be the case. 



APPENDIX A: NOTATION AND CONVENTIONS 

We use units in which h = c = 1, but in which 1, 
so that the Planck mass is given by 



M P = 



1 



32ttG' 



(Al) 



Throughout we use the same sign conventions for metric 
and curvature tensors as in the book of Misner, Thorne, 
and Wheeler ]l8| . Specifically the metric g a b has signa- 
ture (— , +, +, +). Indices k, . . . run over the spatial 
indices 1, 2, 3 while indices a, b, c, . . . run over 0, 1, 2, 3. 
We introduce the metric perturbation 



gab - Vab 



(A2) 



and its trace h = h a t > rj ah , where r\ a b is a flat metric. 
In expressions involving h a b, indices are raised and low- 
ered with the flat spacetime metric rj ab . The coordinate 
derivative of a tensor T a b in a Lorentzian coordinate sys- 
tem with respect to the metric r\ a b is denoted in the usual 
way: 



n c = d c T%. 



The Fourier transform of any 
Minkowski spacetime is defined as 



(A3) 

function F(x) on 



F(k) 



Throughout we use 



x e 



-ikx 



F{x). 



(A4) 



(A5) 



to denote the dot product of two 4-vectors k a and x a in 
Minkowski spacetime. We use Q(x) to denote the step 
function, <d(x) — 1 for x > and otherwise. 
Products of operators A,B are defined as 



(AB)(x,z) = / d 4 yA(x,y)B(y,z). 



(A6) 



Thus, factors of \J—g are not implicit in expressions such 
as GV in our calculations. We show such factors explic- 
itly when they are required, with the exception of the 
notation for products of functions 



{fuh)= / d n x^g(x)h(xyh{x). 



(A7) 
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APPENDIX B: EXPRESSIONS FOR 
CURVATURE INVARIANTS 

In this appendix we expand the various possible local 
counterterms in the effective action to second order in 
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the metric perturbation, and write them in terms of the 
quantity h a b- We define 



h a b=h a b — -^hri abl 
where h — n ab h a b- In n dimensions we then find 



n - 2 



(Bl) 



(B2) 



Note that many of the terms vanish if the Lorentz gauge 
h ab b — is used. In the Lorentz gauge we find that 
(discarding surface terms) 

J d n x8*/=g= J d 4 x(8 -Ah- 2h ab h ab + hh) 

+ ^5 { d A x(Ah - hh) + Q(S 2 ), (BIO) 



and thus 



h a b — h a b — 



n-2 



hrjab- 



(B3) 



In our calculation of the effective action, when we expand 
quantities such as \/—g, R or R ab in terms of h a b and 
h, we find terms of order 0(n — A). Such terms must 
sometimes be kept and not discarded, since they can give 
rise to finite contributions when multiplied by infinite 
terms of the form l/(n — 4) . 

We find for the determinant of the metric tensor 



1 



— !— h - \h cd h cd + 1 
n-2 4 4(n - 2) 



The curvature scalar becomes 

R = R^+R^ + Q(h 3 ), 



where 



rt 1 



1 :h„ a + h; ab 



n — 2 ' 



hh + 0(h 3 ). 

(B4) 

(B5) 
(B6) 



and 



RW 



h h°-b,c ,T Tab c 10 - - 

4 ' A(n-2) 2 

- j- U ~ A l ~h r c h - -h ca ' b h cb a 
(n-2) 2 > c 2 



+ 



' (2hhj ab + h, a h ab l 



n-2\ 
2h ab h 



— h ab h > c 
ac,b " ,b' l ac 



(B7) 



Other useful scalars are 



R 2 = 



1 



(n-2) 

< 2 



h a h b 
2 n ,a a ,b 



n-2 



cd 



h a b ab + o(h 3 ) 



(B8) 



and 



RabR 



ab 



h - c lab d _ V_Z aT b 

n ab c n ,d t n _ 2)2 n ,a n ,b 



2h 



(V 



n — 2 

be j b,ac \ i / )[!,'■> 



h c c n ab + h cab + h cb a 
0(h 3 ). (B9) 



d n x8^R = I d i x(2h ab c c h ab - h/h) 
1 



5 / d /i xh c c h + 0(5 z ), (Bll) 



d n xAj=g~R 2 = I d i xh a a h b b 



and 



6 / <Tx2h a h b " + 0(8'), (B12) 



d n x8^R ab R ab = d 4 x2h ab c c h abd d 



(B13) 



+0(6 2 ), 



where 8 = n — 4. 

Using the Lorentz gauge and specializing to four di- 
mensions, we find for the linearized versions of the local 



curvature tensors ( |4.6| ) and (4.7) 

H ( a l ) (x)^Vabh c c d d -h abc c + 0(h 2 ) 



(B14) 



and 



H a f(x) = -{r, ab h//-h Mc 
+0(h 2 ). 



h c d 
a ab,c d 



(B15) 



APPENDIX C: THE IN-IN PROPAGATOR IN 
MINKOWSKI SPACETIME 



In this Appendix we combine Eqs. (2.31), (p.5[) and 



( 3.8|) s pecial ized to Minko wsk i spacetime to obtain Eqs . 
( 3.13|)-(3.15 ). Using Eq. (3^8) and the expansions (|3.9| ) 
and (3.12), we obtain 



(Fx'A^x^G^x^y) = -S n (x-y)S rt . (CI) 



Note that it follows from the form of A rs (x,y) given by 
Eqs. (EDttl) and (|J) that G°__(x,y) = -G°f + (x,y) and 
that GT^(x,y) = -G°^_(x,y). 
First, the relation 
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G° ++ (x,y) = G(x-y), (C2) 

follows immediately from translational invariance and the 
equation 

- (m 2 - ie)]G° + + (x, y) = -6 n (x - y) (C3) 



which follows from Eq. (CI) 



The equation determining G°_,(x, y) is, from Eq. (CI), 
[a x -(m 2 +ie)}G\(x,y) = 0, (C4) 

which upon Fourier transforming and using translational 
invariance becomes 



-P 



(m 2 + ie)]G°_ + {p) = 0. 



(C5) 



Any function G°_ + (p) which has support only on the hy- 
persurface p 2 = — (m 2 + ie) will be a solution to this 
equatio n, a nd hence G°_ + (p) is not uniquely determined 
by Eq. (C5). What we have not used yet is the boundary 
condition that <f>+ = on the hypersurface given by 
x° = T, where T — > oo. To make use of this boundary 
condition note that the classical equations for cj>± are 

d n yAl s {x,y)<j >s {y) = -J r {x), (C6) 

and that the solutions to these equations are 

4>,{x)= ( d n yG° st (x,y)J t (y)- (C7) 



Enforcing the above mentioned boundary condition 
yields 



d n yG° +s (x,y)J s (y)= d n y G^ s {x,y)J s {y) (C8) 



at x = (T, ), which using Eq. fl2.26| ) simplifies to 

J<Py [G° ++ (x,y) - G°_ + {x,y)] J + (y) 

+ fd n y [G° + _(x,y)-G°__(x,y)] J_(j/)=0. (C9) 

Since the sou rces J±{y) are completely arbitrary, Eqs. 
@ and @ imply that 

G° + (.T,y) = G^-2/) (CIO) 

for a; = (T . . X . X 3 ). Now the Feynman propagator 
G{x) can be written as 

G(x) = -G(x°)A-(x) + 0(-x°)A+(x), (Cll) 

where 



y ' J (2tt)» 



(C12) 

are the pos itive and n egati ve Wightman functions [fl9|| . 
From Eqs. ( |cTc| ) and 0V§ it follows that 

G°_ + (x >y ) = -A-(x-y) (C13) 

for large a; . Equation flC13|) now implies that the ap- 
propriate solution of Eq. flCq) which fulfills the bound- 
ary condition at x a = T corr espon ds to G°_ + (x,y) = 
— A~(x — y), which yields Eq. ( 3.13| ) above. 



APPENDIX D: EXPANSION OF THE 
PROPAGATOR 



In this Appendix we obtain t he e xp ansi on (|3.16| ) for 
the operator G++. From Eqs. (§J), (fH), ( pT2j and 
( |ClD it follows that 

G r t = G° s [S st + V ss 'G s , t + V ss >G s , s „V s »t'G t , t ] 

+0(V 3 ). (Dl) 

This implies that 

y^fO I flO XT /~l0 

++ — U + + "+" U +s V ss'^ s ' + 

+G +s V ss 'G s , s „V s "tG t+ + ... 
= G° ++ [1 + V +S G° S+ + V +s G° ss ,V s , t G° t+ + 0(V 3 )] , 

(D2) 

where we have used the fact that Gj, G5 = 

j4j_,Gj__ = 0. Hence we can write the logarithm of the 
propagator as 



ln(G ++ ) = ln(G° ++ ) + V +S G° S+ + V +s G° ss ,V s , t G° t+ 
-lv +s G s+ V +f G° t+ + O(V 3 ). 



(D3) 



Next we use the fact that V rs is diagonal [cf. Eq. (3.11) 
above] to obtain 

ln(G ++ ) = ln(G° ++ ) + V ++ G° ++ + l V++ G° ++ V ++ G° ++ 
+V ++ G° + _V—G°_ + + 0(V 3 ) 
= ln(G° )+V«G° +^G° 



+Q(V 3 ) , 



^Vg°_^g° 



+ 

(D4) 



which yields Eq. ( 3.1(: ) 



APPENDIX E: TERMS IN THE EFFECTIVE 
ACTION 

In this Appendix we list the terms contributing to 
the effective action ( |3~30| ) given by Eqs. (^32|)-(^39|). 
Throughout we use the Lorentz gauge. We find 



L J 5127T 2 



d*xh(x)(-4m') [Y+- , (El) 



K 2 [h ab 



5127T 2 



d 4 x 



h(x)h(x)(2m 4 ) (Y + 



+2h(x) ab h(x) ab [ -2m 4 Y + - 



(E2) 
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Li[h ab ] = 



1 



5127T 2 

+2h ab {x)h ab (x') 
d 4 k 



d 4 x I d 4 x' ■ 



h(x)h(x') 



(2tt) 4 



Q ik(x—x ) 



46\/ 4 m 2 fc 4 
F+ 15)( m + ~ k +30 



-^-^fc 2 -^(fc 2 +4rn 2 ) 2 ^(fc) 



(E3) 



£2 [h a b] = L 3 [h ab ] 

, . d 4 x / dV<J h{x)h(x') 



V2 / ,4. / ,4,. 



5127T 2 



d 4 fc 



D ik(x—x ) 
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4m 



2m 2 



+ l^! + ^ m 2 (fc 2 +4m2)i?i(fc) 



(E4) 



and 



1 



5127T 2 

d 4 k 



d 4 x / d 4 x'{ h(x)h(x') 



(Y + 4)(2m 4 )-4m 4 F 1 (k) 



ik(x—x') 



Here we have defined the constant 
2 



, 4ir , 
r = +ln— - 7 - In— , 

4 — n /i z /i z 



and the function 
Fi(k) -- 



k 2 + 4m 2 



k 2 



, , arctanh \. 
k 2 \ V k 2 + Am 2 



(E5) 



(E6) 



(E7) 



The term which contains the differences between the 
in-out and in-in effective action is 



U[h a b+,h a b-] 



1 



256tt 2 



d 4 x / d 4 x' 



4m 2 



d 4 k 



ik(x—x ) 



x|/i + (x)/i_(a;')[- 4m 4 + — (fc 2 + 4m 2 ) 

-^(fc 2 +4m 2 ) 2 l Gi(fe) 
+2hf(x)h ab - {x') [ - ^(k 2 + Am 2 ) 2 ] G 1 (fc) I , 



(E8) 



Gi(fc) = 



fc 2 + 4m 2 
fc 2 



9(-fc 2 -4m 2 )e(-fc°). (E9) 
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